Limits of heralded single photon sources based on parametric photon pair generation 
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We derive calculations on the statistics of a heralded single photon source based on parametric 
photon pair generation. These calculations highlight fundamental and practical limits for these 
sources, and show which physical parameters can be optimized to improve the quality of a real 
source. 

PACS numbers: 42.50.Ar; 42.50.Dv 
I. HERALDED SINGLE PHOTON SOURCES II. A SIMPLE HSPS MODEL 



Although the concept of a light particle has existed for 
a long time, the modern description of what would later 
be termed "photon" was introduced by A. Einstein in 
1905 [H. Since then, progress has been made in the de- 
tection of single photons [2| , but the design of true single 
photon sources remains a challenge. One of the proposed 
solutions is a heralded single photon source (HSPS) 1,3 
based on a parametric source of photon pairs. 

We derive herein calculations taking into account rel- 
evant physical parameters of such a source. The simple 
HSPS setup on which these calculations are based en- 
tails no interference effects, so that we can focus on pho- 
ton statistics rather than the full picture of the quantum 
states involved. These simple results shed light onto the 
importance of each individual parameter in the quality 
of the source, and provide formulae that can be used to 
find each parameter's optimum value. They also provide 
frmdamental limits for the performance of HSPS, as these 
sources are inherently imperfect. 

In section [Hi we detail the basic HSPS model. In sec- 
tion |lll] we derive the exact statistics associated with this 
model. In section HVl we examine the properties and lim- 
its of unfiltered HSPS. We then analyze the effect of fil- 
tering in section |Vl We conclude in section IVTl 
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FIG. 1: Simplified HSPS setup. PPPS: parametric photon 
pair source (photons within a single pair are supposed to 
be distinguishable and therefore separable); HB: heralding 
branch; D: detector; EO: electrical output (heralding signal); 
SB: signal branch; OO: optical output (heralded photons). 
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We propose to study the simple setup represented in 
Fig. [TJ In this model, a source of photon pairs (e.g., gen- 
erated via parametric down conversion in a crystal [5| or 
four- wave mixing in a fiber [J] ) is assumed to statistically 
provide N photon pairs with probability P{i-i{N) during 
a specified time bin (e.g., the duration of the pump pulse 
or the active time of a single photon detector) . 

The two photons in each pair are furthermore assumed 
to be distinguishable and separable (e.g., using spatial or 
frequency filtering) . One photon from each pair is used as 
a heralding signal for the second photon. The heralding 
photon is detected by a single photon detector, hence 
announcing the presence of another photon in the signal 
(heralded) branch. For a perfect HSPS, the probability 
of having no photon in the signal branch is reduced to 
zero when a heralding signal is present. 

However, HSPS are fundamentally limited by physics, 
as even a perfect heralding system would not provide a 
true single photon source. Indeed, parametric processes 
inherently generate multiple pairs with a non-zero prob- 
ability. These multiple pairs cannot be eliminated by 
the heralding signal. In fact, we will see below that the 
heralding system increases their probability. In addition, 
none of the physical elements of a real source are perfect. 
In particular, losses occur in both branches and at the de- 
tector level, and dark counts can provide false heralding 
signals. 

Fortunately, it is rather straightforward to model losses 
in a quantum system. Indeed, a lossy system can always 
be modeled by one or several simplified beam splitters for 
which we consider only one input and one output each . 
In addition, some systems (such as the one studied here) 
are so simple that no interference takes place. In those in- 
stances, one needs only to compute intensities (i.e., prob- 
abilities) instead of complex amplitudes for the quantum 
states. This simplification ensures a better readability of 
the results. 

Ill this instance, assuming n-^^ = N photons at the 
input of the simplified beam splitter, < riout £ N 
photons at the output, and a probability of transmission 
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a simple binomial law for the transmission of n photons. 

All losses in a single branch (including detection losses 
due to the limited efficiency of the detector) can be com- 
pounded and modeled by a single beam splitter. We as- 
sume that losses are independent of the mode (i.e., in- 
dependent of the wave-vector direction, spectrum and 
polarization of the photons). Effects of mode-dependent 
filtering are analyzed separately in section IVl 



III. STATISTICS OF THE MODEL 

Let us assume first that N pairs have been produced in 
a single time bin at the input of the system. Let us note 
r]h the transmission efficiency of the heralding branch (in- 
cluding the quantum efficiency of the detector), and rjs 
the transmission efficiency of the signal, or heralded, line. 
Let us also note the probability of a dark count on the 
detector during a single time bin. Using Eq. [U we find 
the probability that a heralding signal is triggered (i.e., 
the detector clicks) to be 
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Remembering that the probability of having N pho- 
tons at the input is noted Pin(iV), the total probabihty 
of having n photons at the output of the signal line, con- 
ditional to the presence of a heralding signal is 
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We can apply Eq.[3]to the most common type of prob- 
ability for the generation of photon pairs, namely Poisson 
statistics. In the case of photon pairs generated via para- 
metric down conversion or four wave mixing, this type of 
statistics arises when a large number of distinguishable 
modes are accessible to the parametric process Q. In 
this case, the A^-photon pair probability is 
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where ^ is the average number of photons pairs generated 
per time bin. 

Injecting Eq. [5] and 2] into Eq. [31 we get 
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where 
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is a correcting term to the Poisson distribution 
6"'^''= (/i?7s)"/n! that would be observed at the output of 
the signal line in the absence of a heralding signal. This 
correcting term usually ensures that the statistics be- 
comes sub-Poisson (as measured by a g^^^ (0) factor lower 
than 1 for the photons at the output of the heralded line). 

There is also a second interesting case where Eq. [3]can 
be applied. When heavy mode filtering takes place down- 
stream of the parametric process, or in the case of some 
inherently narrow processes , the statistics of the pro- 
duced pairs is thermal Q. This case is very interesting, 
because it leads to indistinguishable photons (i.e., pho- 
tons produced in the same spatial and spectral mode), 
which is a requirement for most linear optical quantum 
computing (LOQC) operations Q. The A-photon pair 
thermal probability is 
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Injecting Eq. [2] and [7] into Eq. [3l we get 
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where 
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is a correcting factor to the thermal distribution 
{p.rjsY / (I-|-^77s)"+^ that would be observed at the output 
of the signal line in absence of a heralding signal [ll| . 



IV. PROPERTIES AND LIMITS OF HSPS 
A. Properties 

The form of Eqs. [5] and [6] (rcsp. Eqs. [8] and [H]) makes 
it straightforward to understand the significance of the 
physical parameters of the source. 

First, we have 
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FIG. 2; Example of probability modifications due to herald- 
ing. On the left axis: probabilities associated with the non- 
heralded signal branch (dashed line); improved probabilities 
associated with the heralded signal (full line). On the right 
axis: values of ^p{n) (dash-dotted line). In this example, 
r}h=ris= 50%, dh = 10"" and n = 0.01. 



for thermal statistics respectively. 
For the limit case ^ = we obtain 
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So for 77/j ^ dh (i.e., when the transmission efficiency 
of the heralding branch is much greater than the dark 
count probability on the detector), there is an important 
increase in the probability of seeing one photon versus 
that of seeing no photon in the signal branch. That is 
the main mechanism of HSPS. 

Eqs. |6] and [9] show functions S,p{n) and S,t{n) growing 
with n. This means that multiple pair probabilities in- 
crease even more than the probability of a single pair. 
This is detrimental to the process of getting a true single 
photon source. However, this effect is small compared 
to the attenuation provided by a small ^rjs (see Eqs. [H 
and [8]). Indeed, ^p{n) and ^(ri) arc concave and tend 
asymptotically towards respective finite limits 

lim £n{n) = ^ 



and 



lim 6(") = -T- 

1^+00 (Xh 



In the end, the original exponentially decreasing prob- 
ability (inherent to Poisson and thermal statistics with 
77s M ^ 1) is multiplied by the relatively slowly grow- 
ing correcting factor ^. It is then possible to increase the 
probability for a single photon, while keeping the chances 
of getting multiple pairs low. Fig. [2] shows an example 
of this feature in the case of Poisson statistics. In this 
case, the graphs for thermal statistics would be almost 
identical. 



B. Limits 

Let us examine what happens in a perfect experiment 
where 77s = ?7/i = 1 and dh = 0. 

First, we get ^(0) = 'ft(O) = 0. As expected, all in- 
stances where no photons are present are eliminated. 

Also, the correcting factor for a single photon in the 
signal branch is respectively £,p{l) = e^jie^ — 1) and 
6(1) = (1 + [1)1 [I. Since HSPS are usually operated in 
regimes where 1, these factors can be very large. 

However, the probability of getting exactly one photon 
is not unity, but respectively /i/(e^ — 1) and 1/(1 + /_*). 
These terms become unity only for ^ = 0, i.e., only when 
no photons are generated! This is a fundamental limita- 
tion of HSPS, due to the parametric process used for the 
generation of photon pairs. 

As we have seen, the probability of getting multiple 
pairs is increased by the heralding process. This ef- 
fect can only be managed by reducing [i. However, real 
sources cannot be dimmed too much because of dark 
counts on the detector. 

In order to determine the limits of dimming, we use 
Eq. [5] to compute the photon number variance of a 
Poisson-based HSPS as 
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A true single-photon source would have zero variance. 
In practice, we want the source to exhibit strong sub- 
Poisson behavior, i.e. (A7i)^ ^ (n). 

For fi —>■ 0, dark counts on the detector become dom- 
inant and (Ati)^ (n). The same is true for fi +00, 
this time because the probability of having zero photons 
becomes negligible even in the absence of the heralding 
signal. Hence, minimizing (Ati)^/ (n) provides a value of 
fi that optimizes the sub-Poisson behavior of the source, 
as shown on Fig. [3] for the same set of parameters as 
that of Fig. [2] An equivalent curve for thermal statis- 
tics at the output of the parametric photon pair source 
would exhibit the same behavior, as thermal and Poisson 
statistics are very similar for small values of ^. 



V. EFFECTS OF FILTERING 

A filter (e.g., spectral, spatial or polarization) can be 
placed in one or the other line in order to purify the her- 
alded states (i.e., only one mode is left unfiltered and 
photons in the heralded branch become truly indistin- 
guishable) . 
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FIG. 3: Effect of dimming. (An)'^/ (n) as a function of fi for 
Tih ~ ris = 50% and dh = 10~*. Optimal operation arises for 
~ 0.016 photon per time bin. 



The effect of fihering is different from the effect of 
attenuation. Fihering one of the branches has two ef- 
fects, namely selective attenuation in the filtered branch 
and projection of the photons in the other branch into 
a statistical mixture of photons paired with an unfil- 
tered photon and extraneous photons (formerly paired 
with filtered-out photons). 

In this section, we consider that the pair source fea- 
tures a Poisson statistics, and the filtering selects only 
one mode with a thermal statistics. The extraneous pho- 
tons in the unfiltered branch then retain a Poisson statis- 
tics (i.e., the number of corresponding modes remains 
large). We denote by / the transmitted fraction of pho- 
tons through the filter, so that /i/ is the mean number 
of photons allowed to go through, while /i(l — /) is the 
mean number of photons filtered out of the branch. 



A. Filtering the signal branch 

In order to purify the modes of the heralded photons, 
it is natural to first think about placing a filter in the 
signal branch. In such a configuration, the extraneous 
photons in the heralding branch can simply be seen as 
additional background noise on the detector. It is the 
equivalent of replacing dh by 
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We then have 
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as the filtered mode exhibits thermal statistics. 

The filter improves the modal purity of the source. 
However, comparing equations [8] and [9] with equations fTTI 
and [HI it can be seen that the performance of the source 
is affected by the presence of the filter. In addition to de- 
creasing the number of heralded photons (a feature of fil- 
tering), the filter also increases the noise on the heralding 
line (as dh is replaced by an always larger Vh). This addi- 
tional noise increases the probability of seeing no photon 
in the signal branch when a heralding signal is present, 
which greatly reduces the gain obtained by the heralding 
signal in the first place. 



B. Filtering the heralding branch 

When a filter is placed in the heralding branch, we 
must take into account extraneous photons in the sig- 
nal branch. Hence, the probability of finding exactly n 
photons in the signal branch is 
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is the probability of finding exactly k photons paired to 
an unfiltered heralding photon, and 
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is the independent probability of finding exactly k extra- 
neous photons. 

In the end, we find 
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FIG. 4: Effects of filtering. The original, non-lieraided, prob- 
ability is sliown as a full line. As in Fig. ^ rih ~ rjs — 50%, 
dh = lO""* and fi = 0.01. A filter with / = 10% is then 
inserted on one of the branches. Filtering on the heralding 
branch is shown as a dashed line. Filtering on the signal 
branch is shown as a dash-dotted line. 



which shows that the main contribution to the multi- 
photon statistics is that of the extraneous photons. For- 
tunately, the probability of having multiple photons in 
the signal branch remains negligible for sufficiently small 

Fig. m shows how filtering in the heralding and signal 
branch modifies the probabilities of the source. All filter- 
ing in a single branch is both beneficial to the purity of 
the modes and detrimental to the sub-Poisson quality of 
the source. The physical effects of filtering in the signal 
and heralding branches are different. Filtering in the sig- 
nal branch degrades the source by increasing Ps(0), which 
goes against the heralding mechanism of the HSPS. On 
the other hand, a filter placed on the heralding branch 
increases multiple pair probabilities. This effect can be 
reduced by proper dimming of the photon pair source. 
Hence, filtering the heralding branch is a better solution 
for HSPS. 
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and Cn is the n*^ order Laguerre polynomial defined 
as [H Eq. 22.11.6] 



fe=0 



k k\ 



(26) 



As expected, Eq. [23] reduces to Eq. [9]for / = 1. How- 
ever, for / < 1, ^/t(n) no longer features a finite limit 
when n —>■ +oo {^h grows exponentially), and multi-pair 
probabilities increase. Indeed, as n grows, Eq. [22] be- 
comes 
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VI. CONCLUSION 

We have developed a model for heralded single pho- 
ton sources based on parametric photon pair generation. 
This model shows that the heralding process results in a 
multiplying factor that increases the probability of seeing 
at least one photon at the optical output when a herald- 
ing signal is present. 

Because of the parametric photon pair source used in 
HSPS, the probability of seeing multiple photons at the 
output cannot be eliminated. In fact, the heralding pro- 
cess only increases that probability. Hence, the probabil- 
ity of having exactly one photon is never unity. The only 
way to reduce the occurrence of multiple pairs is to lower 
the average number of pairs produced by the parametric 
source. However, dark counts on the detector impose a 
lower limit on that average value. 

When a modally pure source of photons is required, 
the best solution is to try and increase the purity of the 
parametric photon source. However, this is not always 
possible, and additional filtering might be required inside 
the heralding system. In this case, it is preferable to 
implement filtering of the heralding line, rather than the 
signal line. 
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